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We present a quantitative estimate of the anisotropic power and scaling of magnetic field fluctu- 
ations in inertial range magnetohydrodynamic turbulence, using a novel wavelet technique applied 
to spacecraft measurements in the solar wind. We show for the first time that, when the local 
magnetic field direction is parallel to the flow, the spacecraft-frame spectrum has a spectral index 
near 2. This can be interpreted as the signature of a population of fluctuations in field-parallel 
wavenumbers with a k7 spectrum but is also consistent with the presence of a "critical balance" 
style turbulent cascade. We also find, in common with previous studies, that most of the power is 
contained in wavevectors at large angles to the local magnetic field and that this component of the 
turbulence has a spectral index of 5/3. 

PACS numbers: 52.35.Ra,96.50.Bh,52.30.Cv,95.75.Wx 



Magnetised plasmas fill most of the Universe and in 
many regions, turbulence plays important roles in the 
transport of energy and momentum and the accelera- 
tion and scattering of charged particles. Many aspects of 
plasma turbulence remain poorly understood, however. 
Here we present results on one of these, the anisotropy 
of the energy spectrum of magnetohydrodynamic (MHD) 
turbulence with respect to the magnetic field. 

In classical hydrodynamics, velocity fluctuations Suk 
with a wavenumber k decay and transfer energy to 
smaller scales on the shear timescale, rg ^ l/{kSuk). 
Within the steady inertial range, far from the energy in- 
put ("outer") and dissipation scales, this leads to the 
dimensional result (Suk)^ oc e/fc, where e is the energy 
dissipation rate per unit mass. This gives the familiar 
Kolmogorov energy spectrum P{k) oc k~^/^, widely ob- 
served in hydrodynamic turbulence. In a plasma, fluc- 
tuations can also propagate, as Alfven waves parallel to 
the magnetic field, and this leads to the Alfven timescale, 
TA ~ l/(fcj|V}i), being dynamically important. Here fcy is 
the component of the wavevector of the fluctuation par- 
allel to the local magnetic field and Va the Alfven speed. 
If TA <S Ts and assuming isotropy with respect to the 
local field, this leads to Iroshnikov-Kraichnan turbulence 
[l|,i] where (SukY ex eVk//c and P(fc) oc k~^^^ [e.g.ii]. 

However, measurements in both space plasmas and ter- 
restrial plasma devices have shown that turbulent fluc- 
tuations are not isotropic. They typically have much 
longer correlation lengths along the field than across it 



[a 0) S S M, [i3l and the spectral index for the mag- 
netic energy is nearer 5/3 than 3/2 [111, Il4|- When 
there is an energetically significant large-scale magnetic 
field, anisotropic models of MHD turbulence are required 
[l^. M. [Tsl [Ta. [TtI [ill . For example, in the "critical bal- 
ance" framework [14| . turbulent energy evolves towards 
wavevectors where the shear and Alfven timescales are 
balanced and most power resides in wavevectors where 
TS < TA, i.e. fc|| < k]^\^/^V^\ 

The solar wind is a unique environment in which to 
study space plasma turbulence: it is relatively accessible 
and can be directly measured in exquisite detail using 
spacecraft instruments [e.g. [ij, \1m 0, |2i|. The solar 
wind flows radially away from the Sun at a velocity V 
of several hundred kms~^, much faster than spacecraft 
motions (a few kms~^) or the plasma wave speeds (tens 
of kms~^). As a result, in the plasma frame spacecraft 
measure along a radial line. Using Taylor's hypothesis 
|2a | , one can relate the spacecraft frame energy spectrum 
■p(/) to the wavevector spectrum P{k) |23j : 



^(/) = / d^k P{k) 6{2nf -k-V). 



(1) 
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Anisotropics in P{k) with respect to the magnetic field 
can be analysed by measuring how V{f) varies with the 
angle of the magnetic field to the flow, 6b- 

The exact form of this anisotropy is unknown, but one 
can make approximations motivated by theory and com- 
pare predictions with observations. One simple approx- 
imation is to assume that P{k) — except for wavevec- 
tors exactly parallel (so-called "slab") or perpendicular 
( "2D" ) to the local magnetic field [2J| . The correspond- 



ing frequency spectrum can be deduced from Eq. [TJ 
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where Cgiab and C2D sue constants and 'jsiab and 72D are 
the spectral indexes of these components. V is thus in- 
sensitive to the slab component when the field is perpen- 
dicular to the flow and insensitive to 2D when the field 
is parallel [2J|. Crucially, one can determine the scaling 
of both components by measuring the spectral index a 
oiV Pi /-" separately for Ob = 0° and Ob = 90°. 

In the case of critical balance [ij] , the 3D power spec- 
trum takes the form 



P{k) oc kl^^'^g 



VAk\\ 



£l/3fc 



2/3 



(3) 



V{f] Ob) then depends in a complicated way on the un- 
specified function g{y), but one can show that Eq. [3]im- 
plies P{f;eB = 90°) « /"^/s and V{f;0B = 0°) « f-^ 
with the latter also smaller in magnitude at a given /. 
This result is independent of the precise form of g{y). 
In the case of critical balance we would therefore expect 
an anisotropy in both the power levels and the spectral 
index of the spectrum. 

Here, we use 30 days (1995, days 100-130) of 1 second 
resolution measurements of magnetic field fluctuations 
[25| by the Ulysses spacecraft. During this time, Ulysses 
was within the steady high speed {V ~ 750kms~^) so- 
lar wind from the Sun's Northern polar coronal hole at 
1.4 AU from the Sun [20|. Fluctuations within the so- 
lar wind inertial range, corres pon ding to spacecraft time 
scales of seconds to minutes [27[, are superimposed on 
large amplitude ( | (5.B | / 1 i? | ~ 1) Alfven waves on time 
scales of hours [28| which result in large variations in 
6b ■ The minimum variance direction of the inertial range 
fluctuations follows the local magnetic field direction very 
closely [29,], indicating that the local field orders the 
behaviour of the fluctuations. We can therefore study 
the anisotropics of the turbulence by measuring how the 
spacecraft frame spectrum of magnetic fluctuations varies 
with 9b- We perform this analysis using a new wavelet 
method, which is sensitive to the constantly changing lo- 
cal magnetic field direction. 

Wavelets have been used extensively to study physical 
time series [e.g. |30, [Ml IH, ^M- The Morlet wavelet is 
relatively well-localised in frequency, being rather wave- 
like fH: 
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With tiJo = 6, it is possible to construct a nearly orthonor- 
mal set of wavelets. For each magnetic field component 
i, we have a time series Bi(tk), where ife = tg + kSt and 
6t = \s. The discrete wavelet transform wi (tj , /;) of such 
a time series, at a time tj and frequency // is given by 
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The time scale or dilation parameter si is directly re- 
lated to the peak frequency response /; of the wavelet. 
For ujQ = 6^ Sl = 1.031/// [e.g. [32]. In practice, the 
wavelet transform is more efficiently calculated using 
Fourier transforms rather than directly in the time do- 
main. We calculate the wavelet coefficients at ten log- 
arithmically spaced frequencies fi = fo ■ (8/5)"', where 
/ = 0, 1, . . . , 9: /o = 0.25 Hz and /g = 3.6 mHz. 

The wavelet coefficients Wi can be used to calculate the 
power in a time series: at a time tj and frequency // in 
component i the power is proportional to 



Vu{t,Ji)^fi\w,it,Ji)\' 



(6) 



Here we analyse the trace, V = 'SVu- 

Fluctuations at a given scale are sensitive to the local 
magnetic field - but the definition of "local" varies with 
the spatial scale of the fluctuations of interest. In general, 
one would expect fluctuations with a given wavelength to 
be sensitive to the magnetic field on approximately this 
scale and above. 

In order to measure the scale-dependent local magnetic 
field direction, we calculate the amplitude envelope of the 
Morlet wavelet, \tp{tj, fi)\'^ - this is a Gaussian centered 
on time tj with a width of 1.67s; - and calculate the sum 
over the data set of the product of this envelope with the 
magnetic field time series, for each field component i: 
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For each s/, this results in a time series of vectors b{tj,si) 
which point in the direction of the local magnetic field 
at the time tj , associated with a time scale si (or equiva- 
lently /;). This provides a frequency- and time- localised 
mean field direction, and hence 9b and the azimuthal 
angle 0s, for every wavelet coefficient. 

We next construct a set of 404 bins, each of which 
subtends approximately the same solid angle, and which 
together cover all directions. For each bin, for a given fre- 
quency /; and field component i, we select all wavelet co- 
efficients that have average magnetic field angles {9b,<Pb) 
within the bin. The mean of these coefficients is then the 
average power in component i at frequency // when the 
field points in that direction. 

This process results in a scale-sensitive estimate of the 
magnetic field power spectrum as a function of the mag- 
netic field angle relative to the solar wind flow (sampling) 
direction - that is, we estimate 'P(/; 0_b, c/jb)- Many bins 
contain thousands of measurements; in order to ensure 
reliable statistics, any angle bin with fewer than 40 con- 
tributing power levels is rejected. If the fluctuations are 
axisymmetric around B, V should be independent of 4>b- 
Our measurements indicate that this is indeed the case 
and we therefore consider values averaged over all 4>b- 

Using this wavelet method, we can estimate the 
spacecraft-frame power spectrum at a range of frequen- 
cies / and magnetic field/flow angles 9b ■ Fig. [T] shows 
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FIG. 1: Magnetic power spectra at two different angle ranges 
of the local magnetic field to the flow: 0-10° (circles) and 80- 
90° (diamonds). Note the reduced power levels and steeper 
slope associated with the smaller angle. Guide lines with 
slopes of 5/3 and 2 are shown above and below the data. 
Spectral indices in Fig. [5] are calculated over the scales be- 
tween the dotted vertical lines. 



typical power levels for two ranges of Ob'- 0°-10° and 80°- 
90°. Although both are well described by power laws in 
/ over the range of frequencies considered here, it is ap- 
parent that the power levels for 6b = 0°-10° are lower 
than those for 80°-90°, in agreement with several previ- 
ous studies [2J, |3J, Is^ and consistent with the expecta- 
tion that most power in the fluctuations is in wavevectors 
at large angles to the magnetic field [l\. In addition, the 
spectrum for 6b = 0°-10° is steeper than that for larger 
angles. 

The variations in the spectrum power level and spectral 
index with 6b are more easily seen in Fig. [5J where it 
is apparent that there is a smooth variation in power 
with field/flow angle. Note that occasional folds in the 
magnetic field past 6b = 90° (due to the presence of large 
amplitude Alfven waves [25|) mean that it is possible to 
measure variations for 6b > 90°, although not all angles 
can be measured. 

The most important result in this paper is shown in 
the bottom panel of Fig. ^ where it is clear that there 
is a systematic variation in spectral index with 6b- For 
most angles, a « 5/3, in accord with most previous solar 
wind measurements [e.g.ll2|. However, when 0b — > 0°, a 
changes in a smooth manner towards a value of around 
2. A spectral index of 2 at small field/flow angles has 
not previously been reported. It is strong evidence of 
anisotropic energy transfer in the MHD cascade and 
was predicted for any "critical balance" type of cascade 
[la . ISq . A spectral index of 2 at 6b = 0° and 5/3 
at 9b = 90° is consistent with the presence of a criti- 
cal balance cascade. However, it is surprisingly difficult 
to distinguish between the critical balance and slab/2D 
approximations from these results and indeed they are 
also broadly consistent with a dominant population of 2D 
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FIG. 2: Top panel: Trace of power in the magnetic fleld as a 
function of the angle between the local magnetic field and the 
sampling direction at a spacecraft frequency of 61 mHz. The 
larger scatter for 6b > 90° is the result of fewer data points 
at these angles. Bottom panel: spectral index of the trace, 
fitted over spacecraft frequencies from 15-98 mHz. 



fluctuations with 72D — 5/3 and a smaller slab popula- 
tion with "fsiab — 2. The small range of scales over which 
we measure the fluctuations does not rule out other scal- 
ings (e.g. exponential) for a possible low amplitude slab 
component. 

It is perhaps surprising that the f~^ scaling has not 
previously been observed in the solar wind and indeed we 
are aware of ongoing work by others which does not show 
such variation. However, the constantly changing back- 
ground field direction means that only by using very short 
averaging periods can we avoid "smearing out" the very 
low power fluctuations observed when 6b ~ 0, which is 
the only time that we see the f~^ scaling. Recent multi- 
spacecraft studies (which are not susceptible to these ef- 
fects) have not revealed this scaling in the solar wind 
[37| . although anisotropic scaling has been observed in 
the magnetosheath |38| . 

Note that the variation of P{6b) in Fig. [5] is not sym- 
metric around dB = 90°, which is not possible under 
our assumptions of homogeneity and Taylor's hypothe- 
sis, and therefore one or both of these assumptions must 
to some extent be violated. The lack of symmetry in 
Fig. [5] may be due to kinetic effects at the small scale 
edge of the inertial range, although Fig. [T] suggests that 
the steeper spectrum for small 6b extends over a wide 
range of frequencies. The dynamical effects which cause 
local changes in the field direction, including large-scale 
Alfven waves and microstreams, might cause systematic 
changes in power levels with 6b- However, while they 
might change P{6b), they would be unlikely to change 
the a: we consider the measurement of a steeper spec- 
trum near 6b = 0° to be a robust result. 



We are only measuring the magnetic energy spectrum 
here, and not that of the plasma velocity: telemetry lim- 
itations make it impossible to study the velocity fluc- 
tuations on these timescales using Ulysses data. How- 
ever, other observations [ll| indicate that the solar wind 
magnetic and kinetic energy spectra often have distinct 
slopes, typically 5/3 and 3/2 respectively. Many sim- 
ulation studies also find distinct slopes, as does a re- 
cent closure model [33 , but most theoretical models give 
forms for the total (kinetic plus magnetic) energy spec- 
trum, making it difficult to compare them directly with 
our results. Nonetheless, our results are in accord with 
simulations where the amplitude of the background field 
is approximately the same as the rms B- as occurs in 
high-latitude solar wind [40| - and which yield kinetic 
and magnetic spectra with slopes of ~ 5/3. Our results 
are not consistent with spectral indexes of 3/2 associ- 
ated with strong background field simulations and mod- 
els[e.g.,[l5, 39, 41]. 

The anisotropy reported here may also be influenced 
by the anisotropic "driving" of the turbulence at large 
scales by large-amplitude, predominantly anti-Sunward- 
propagating, Alfven waves. Such driving is in contrast 



to the probably isotropic and possibly weak-amplitude 
injection occurring in many astrophysical plasmas, and 
may limit the applicability of our findings to other plasma 
regimes. 

The range of 6b over which the spectral index deviates 
from 5/3 is rather larger than we would expect on the 
basis of the variation of V with 9b , under either a critical 
balance or slab/2D cascade framework. Currently we 
do not have a good explanation for this discrepancy but 
hope to address it in a later paper. 

Finally, we note that the wavelet method can be ex- 
tended beyond just the trace of the spectrum as used 
here. We can also measure the individual elements of the 
power spectral tensor, revealing information about the 
field-parallel variance (diagonal elements), helicity (off- 
diagonal elements) and even intermittency (using higher- 
order moments). We intend to present our analyses of 
these measurements in the near future. 

Ulysses analysis at Imperial College London is sup- 
ported by the UK STFC. The authors are grateful to 
ISSI, Bern, for their support of this work. T. Horbury ac- 
knowledges helpful conversations with A. Schekochihin, 
K. Osman and C. Chen. 
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